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Abstract 

The recent measurements of the color-suppressed modes — > D^*'^^7r^ unply non-vanishing 
relative final-state interaction (FSI) phases among various B Dtt decay amplitudes. Depending 
on whether or not FSIs are implemented in the topological quark-diagram amplitudes, two solutions 
for the parameters ai and 02 are extracted from data using various form-factor models. It is found 
that |a2(-D7r)| 0.35 — 0.60 and |a2(-D*vr)| ~ 0.25 — 0.50 with a relative phase of order 60° between 
fli and 02. If FSIs are not included in quark-diagram amplitudes from the outset, /af^ and 

will become smaller. The large value of |a2(D7r)| compared to \a^{DTT)\ or naive expectation 
implies the importance of long-distance FSI contributions to color-suppressed internal VF-emission 
via final-state rescatterings of the color-allowed tree amplitude. 



PACS numbers: 13.25.-k 
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I. INTRODUCTION 



For some time B J/ipK and B J/ipK* remain to be the only color-suppressed B 
meson two-body decay modes that have been measured experimentally. Recently, the long 
awaited color-suppressed decay modes B^ Z}(*)°7r° are finally measured by both Belle 
and CLEO [@] with the D^vr^ branching ratio larger than the upper limit previously reported 
0]. The channels B° ^ D^*^^t] and B° D^*^°u are also observed by Belle [|]. We shall 
see below that the theoretical predictions based on the factorization approach in general are 
too small to account for the observed decay rates of color-suppressed modes D^*^^X'^ with 
X = 7c,ri,uj. This has important implications for final-state interactions (FSIs). 

Under the factorization hypothesis, the nonleptonic decay amplitudes are approxi- 
mated by the factorized hadronic matrix elements multiplied by some universal, process- 
independent effective coefficients Oj. Based on the factorization assumption, one can catalog 
the decay processes into three classes. For class-I decays, the decay amphtudes, dominated by 
the color-allowed external l^-emission, are proportional to ai(Oi)fact where Oi is a charged 
current-charged current 4-quark operator. For class-II decays, the decay amplitudes, gov- 
erned by the color-suppressed internal l^-emission, are described by a2(02)fact with O2 being 
a neutral current-neutral current 4-quark operator. The decay amplitudes of the class-Ill 
decays involve a linear combination of ai(Oi)fact and a2(02)fact- If factorization works, the 
effective coefficients Oj in nonleptonic B or D decays should be channel by channel indepen- 
dent. 

What is the relation between the coefficients and the Wilson coefficients in the effective 
Hamiltonian approach ? Under the naive factorization hypothesis, one has 

= Cl(/i) + ^C2(/i), a2(/^) = C2(^) + ^Ci(/i), (1) 

for decay amplitudes induced by current- current operators Oi,2(At), where ci,2(At) are the 
corresponding Wilson coefficients and Nc is the number of colors. In the absence of QCD 
corrections, Ci = 1 and C2 = 0, and hence class-II modes governed by 02 = l/Ai'c are 
obviously "color-suppressed". However, this naive factorization approach encounters two 
principal difficulties: (i) the coefficients given by Eq. (0) are renormalization scale and 
75-scheme dependent, and (ii) it fails to describe the color-suppressed class-II decay modes. 
For example, the ratio R = r(D° 'K^tt^)/T{D'^ K^-n^) is predicted to be only of order 
3 X 10^'^ due to the smallness of 02 in the naive factorization approach, while experimentally 
it is measured to be 0.55±0.06 It is known that the decay —>■ K*7r° is enhanced by two 
mechanisms. First, 02 receives a large nonfactorizable correction. Second, the weak decay 
— > K^n^ followed by the inelastic rescattering K^n^ K^t:'^ can raise i3(Z^° — > K^7r°) 
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dramatically by lowering B{D^ K vr"*"). 

Beyond naive factorization the parameters ai.2 have the general expression 

^1,2 = C2,i(yu) H — ^^77^ + Honfactorizable corrections, (2) 

where nonfactorizable corrections include vertex corrections, hard spectator interactions in- 
volving the spectator quark of the heavy meson, and FSI effects from inelastic rescattering, 
resonance effects, ■ ■ ■, etc. In the generalized factorization approach of one includes the 
vertex corrections which will compensate the renormalization scale and 75-scheme depen- 
dence of the Wilson coefficients to render ai 2 scale and scheme independent. Contrary to 
the naive one, the improved generalized factorization scheme assumes that nonfactorizable 
effects are incorporated in a process independent form. Since not all nonfactorizable effects 
are calculable by perturbative QCD, one will treat ai and 02 as free parameters in the gen- 
eralized factorization approach and extract them from experiment. The phenomenological 
analysis of two-body decay data of D and B mesons will tell us if the generalized factoriza- 
tion hypothesis works reasonably well by studying the variation of the parameters 01^2 from 
channel to channel. 

The experimental measurement of 5 ^ J/ipK leads to \a2{J /ipK)\ = 0.26 ± 0.02 [0. 
This seems to be also supported by the study of 5 ^ Dn decays: Assuming no relative 
phase between ai and a2, the result 02 ~ (9(0.20 — 0.30) ||^J^ is inferred from the data of 

^ D^*)+n- and D^*^^n-. However, as we shall show below, the above value of 

02 leads to too small decay rates for 5^ D^^'^^n^ when compared to recent measurements. 
In order to account for the observation, one needs a larger a2{D'K) with a non-trivial phase 
relative to Oi. The importance of FSIs has long been realized in charm decay since some 
resonances are known to exist at energies close to the mass of the charmed meson. We shall 
see in this work that, just as K^tt^, both nonfactorizable effects and FSIs are also 

needed to explain the data of 5^ — > D^*^^^^, though these two effects in B decays are naively 
expected to be not as dramatic as in the charm case. 

The color-suppressed mode is a very suitable place for studying the effect of FSIs (es- 
pecially the soft one) in weak decays. The ratio of the color-suppressed decay amplitudes 
with and without FSIs is %^ = \A{D^ K^n^)/A{D^ ^°7r°)^ithout fsisI ~ 2.0 and the 
relative phase between — > K^n'^ and K^n^ is about 150°. It is expected that for 

B DiT decay, Rd-k and the relative phase among decay amplitudes will become smaller. 
The recent measurement of the B^ D^n^ mode allows us to determine the above two 
quantities. We shall see that although the relative phase among B Dn decay amplitudes 
becomes smaller, -R^ic.)^ does not decrease in a significant way from charm to bottom case. 
The implications and related physics will be discussed below in details. 
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II. FACTORIZATION 



We begin with by considering the branching ratios of the color-suppressed modes B 
]j{*)^X^ {X = TT, 7], (jj) within the framework of the factorization approach. The ^ D^n^ 
amphtude is given by 

^(#^DV) = i=(-C + £), (3) 

where C, £ are color-suppressed internal VF-emission and VF-exchange amplitudes, respec- 
tively. In terms of the factorized hadronic matrix elements, they read 

^c6K:.«2(^7r)(m^ - m^)/BF°-^-(m|), (4) 

where a2{DTT) is a parameter to be determined from experiment. The annihilation form 
factor FQ^^'^{m\) is expected to be suppressed at large momentum transfer, = m^, 
corresponding to the conventional helicity suppression. Based on the argument of helicity and 
color suppression, one may therefore neglect short-distance (hard) W^-exchange contributions. 
However, it is not clear if the long-distance contribution to l^-exchange is also negligible. 
Likewise, 

D'tj) = V^,V:,a,{Dr^){ml - ml)fnF^\mll 
^ DV) = -^Vc,V:,a2{D*7r)V2mn*fn*Ff-{ml.), (5) 



and 



rriB + ^ tub + 



(6) 



Here factorization implies a universal 02, namely, a2{D*uj) — a2{Duj) — a2{Dr]) — a2{D*Tr) — 
a2{D7r). In naive factorization, 02 is not only small, of order 0.10, but also renormaliza- 
tion scale and scheme dependent. In the generalized factorization approach, the scale- and 
scheme-independent 02 can be extracted from experiment and the factorization hypothesis 
is tested by studying 02 to see if it is process independent or insensitive. 
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To proceed, we shall consider four distinct form-factor models: the Neubert-Rieckert- 
Stech-Xu (NRSX) model 0, the relativistic light-front (LF) quark model [0], the Neubert- 
Stech (NS) model [§, and the Melikhov-Stech (MS) model based on the constituent quark 
picture Q. The NRSX model takes the Bauer-Stech-Wirbel (BSW) model |12| results 
for the form factors at zero momentum transfer but makes a different ansatz for their 
dependence, namely, a dipole behavior is assumed for the form factors Fi, Aq, A2, V, 
motivated by heavy quark symmetry, and a monopole dependence for Fq,Ai, where we have 
followed the definition of form factors given in [jl2|. For reader's convenience, the values of 
relevant form factors are listed in Table I (see 0] for some details about the NS model). 

The form factors for B —>■ 7] and B —>■ r]' transitions have been calculated by BSW [jl2l in 
a relativistic quark model. However, in their relativistic quark model calculation of B rj^'^ 
transitions, BSW considered only the uu component of the rj and rj'; that is, the form factors 
calculated by BSW are actually Fq^"""" and Fq^""^ induced from the h ^ u transition. It is 
thus more natural to consider the flavor basis of rjq and rjs defined by 

r]q = -^(mm + dd), rjs = ss. (7) 

The wave functions of the 77 and rj' are given by 





where (p = 9 + arctanv^, and 9 is the 77 — 77' mixing angle in the octet-singlet basis. The 
physical form factors then have the simple expressions: 



=^^os^ F,^r , <r = ^ sin F^<^ . (9) 



Using F(f''""(0) = 0.307 and Fo^''""(0) = 0.254 obtained from |T2[ and the mixing angle 
= 39.3° (or 9 = -15.4°) we find Fq^^^O) = 0.168 and Fo^'''(0) = 0.114 in the BSW 
model and hence the NRSX model. For other form-factor models,^ we shall apply the relation 
based on isospin-quartet symmetry 

pBv^u ^ fI^P'-^ = F^J (10) 
and Eq. (^ to obtain the physical B — rj and B — rj' transition form factors. 



*The form factors FQ^{m'\)) = 0.28 and F^^{nn?j^,) = 0.33 for the NS model obtained in are 
larger than ours by about a factor of 2. 
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TABLE I. Form factors in various form-factor models. Except for the NRSX model, the rela- 
tions ^f'^(g^) = Af^\q'^) {i = 0, 1,2) and V^'^{q'^) = V'^p\q^) are assumed in all the form-factor 
models. The pion in the B — tt transition is referred to the charged one. 







) F^^ml, 


) n^"!™?, 


) Ff^iml, 


) Fi'^irn 




) Af-iml, 


) Ai-{ml, 


) V^-iml,) 


NRSX 


0.37 


0.45 


0.19 


0.23 


0.69 


0.62 


0.26 


0.23 


0.27 


0.32 


LF 


0.34 


0.39 


0.18 


0.22 


0.70 


0.73 


0.25 


0.16 


0.17 


0.28 


MS 


0.32 


0.36 


0.17 


0.20 


0.67 


0.69 


0.26 


0.20 


0.21 


0.28 


NS 


0.27 


0.32 


0.15 


0.18 


0.63 


0.64 


0.22 


0.20 


0.22 


0.22 



As mentioned in the Introduction, in the absence of a relative phase between ai and a2, 
a value of 02 in the range of 0.20 to 0.30 is inferred from the data of — > D'^*^+7r^ and 
D^*^^^^ . For definiteness, we shall use the representative value 02 = 0.25 for the 
purpose of illustration. The calculated branching ratios for 5^ d(*)ox^ are shown in 
Table II for ffj = 200 MeV and fo* = 230 MeV. Evidently, the predicted rates for color 
suppressed modes are too small compared to recent measurements. It should be stressed 
that if there is no relative phase between ai and 02, then one cannot increase 02 arbitrarily 
to fit the data as this will enhance the decay rate of the A/ = 3/2 mode B^ D^*^^n^ and 
destroy the agreement between theory and experiment for the charged mode. For example, 
fitting 02 to the data of D'^tt^ without FSIs will yield 02 = 0.45 in the MS model, which in 
turn implies B{B^ D^n^) = 7.9 x 10^^ and this is obviously too large compared to the 
experimental value (5.3 ± 0.5) x 10^'^ [Q. In this case, one needs FSIs to convert D^n^ into 
D^n^. In contrast, if 02 is of order 0.45, then a relative strong phase between ai and 02 will 
be needed in order not to over-estimate the D^7t~ rate. In either case, we conclude that FSIs 
are the necessary ingredients for understanding the data. 



III. EXTRACTION OF m AND 02 



In this section we will extract the parameters ai and 02 in two different approaches. In 
the first approach, the topological amplitudes are assumed to incorporate all the information 
of strong interactions. Therefore, ai_2 thus determined already include the effects of FSIs. In 
the second approach, one will assume that quark-diagram topologies in their original forms 
do not include FSIs from the outset. 
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TABLE II. Predicted branching ratios (in units of 10"^) of d(*)Ox° {X = 7r,ri,uj) in the 

generahzed approach with various form-factor models for 02 = 0.25, fo = 200 MeV and Jd* = 230 
MeV. 



Decay mode 


NRSX 




1.13 


B° ^ Z)*o^o 


1.57 




0.55 




0.76 




0.76 




1.60 



LF MS NS 



Experiments 



Belle M CLEO |2| 



0.93 0.82 0.58 3.1 ±0.4 ±0.5 2.74l°J^ ± 0.55 

1.20 1.01 0.80 2.7l[]:^t|]| 2.20l|];i ± 0.79 

0.53 0.48 0.34 1.4+[j:| ± 0.3 

0.68 0.58 0.46 2.0l[5| ± 0.4 

0.71 0.76 0.54 1.8±0.5l[5j 

1.16 1.75 1.35 3.ll{:?±0.8 



A. Direct analysis 

In terms of the quark- diagram topologies T, C and S, where T is the color-allowed 
external ly-emission amplitude, the other B Dtt amplitudes can be expressed as 

A(B^ D+71-) = T + S, 

A{B- ^ D\-) =T + C, (11) 

and they satisfy the isospin triangle relation 

A(B^ D+TT-) = V2A(^ DV) ± A{B- D\-). (12) 

In writing Eqs. (^ and (|l5) it has been assumed that the topologies T, C, £ include 
the information of all strong interactions for physical B Dtc amplitudes (for an earlier 
discussion of quark-diagram amplitudes, see ||15|). Now since all three sides of the B 



Dtt triangle are measured, we are able to determine the relative phases among the decay 
amplitudes. Using the data 0] 

i3(5° ^ D+Ti-) = (3.0 ± 0.4) X 10-^ B{B- D\-) = (5.3 ± 0.5) x 10'^ 
i3(#' ^ D*+7r") = (2.76 ± 0.21) X 10"^ -B(5" ^ D*°7r") = (4.6 ± 0.4) x 10~^ (13) 

and the combined value of Belle and CLEO for the neutral modes (see Table II) 

B(B° D°7r°) = (2.92 ± 0.46) x 10-^ B(B° D*%°) = (2.47 ± 0.67) x 10-^ (14) 
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we find (only the central values for phase angles are shown here) 



C-£ 



C-S 



Dn 



r + £ 



(0.44 ± 0.05) e 
(0.42 ± 0.06) e 



C-S 



i63° 



r + c 

C-S 



DiT 



(0.34 ±0.03) e 



i37° 



(0.34 ±0.05) e 



i44° 



(15) 



D*n T+C 

where we have employed the B meson lifetimes given in 0]. 

The same phases also can be obtained from the isospin analysis. Decomposing the phys- 
ical amplitudes into their isospin amplitudes yields 



A{B- - 



D- 



TT 



■ DOtt") 
D\-) 



Al/2 + \l ^^3/2, 
^3/2, 



l3/2- 



V 3 

V3A3/2- 



'1 

l2 



(16) 



The isospin amplitudes are related to the topological quark-diagram amplitudes via 



A,/2 = ^{2r-C + 3S), 



A 



3/2 



(17) 



Intuitively, the phase shift difference between A1/2 and A3/2, which is of order 90° for D 
Ktt modes (see below), is expected to play a minor role in the energetic B Dtt decay, 
the counterpart of D ^ Kit in the B system, as the decay particles are moving fast, not 
allowing adequate time for final-state interactions. Applying the relations (see e.g. @]) 

1 



\A 



1/2 1 



\AiB^ 



D+n-)\^ + \A{B 



0^0m2 



\A 



3/2 



cos 



{5i/2 — ^3/2) 



'i\A{B^ 



r = ^iA(i?" 



.D°vr-)p, 



2|Ai/2 



\A 



3/2 1 



2V2IA1/2IIA 



I3/2I 



we obtain 



Al/2 



72^3/2 



(0.70 ± 0.10) e 



i29° 



A, 



/2 



Dtt 



V2A3/2 



(0.74 ± 0.07) e 



i29° 



(19) 



D*TT 



Similar results are also obtained before in []TB|,|T7| using the preliminary Belle and CLEO 
measurements. It is easy to check that the ratio {C — S)/{C + S) in Eq. ( [T^ ) follows 
from Eqs. ( p!7D and (|19|). It is also interesting to compare the above results with that for 
D^K^*\ decays §: 

- (2.70 ± 0.14) e^^°° 



V2A 



3/2 



Ktt 



V2A 



3/2 



(3.97 ± 0.25) e 



jl04° 



(20) 



K TT 
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The smaller isospin phase shift difference in B decays is in accord with expectation. Notice 
that while A/ = 1/2 and 3/2 amplitudes in B ^ D^*^7r are of the same size, the D Kn 
decays are dominated by the isospin AI =1/2 amplitude. In the heavy quark limit, the 
ratio of Ai/2/{V2A^/2) approaches to unity [T^. Evidently, the charm system exhibits a 
more deviation than the B system from the heavy quark limit, as expected. 

The ratio of 02/01 can be extracted from Eq. (p^Sf) or Eq. (plOl). Noting that the factorized 
color-allowed tree amplitude reads 



T 



I— VcbV*^ai{D'K){m 



V2 



(21) 



and neglecting VT-exchange contributions, we get 



02 

Oi 



(0.44 ± 0.05) e 



j59° 



X 



m 



D 



?BDi 



ml) 



DiT 



fo ml - ml F^^'iml) 



1 - 


Al/2 






D-K 


\+ 








Dn 



fD ml - ml F^^[ 



{22) 



m 



D) 



Likewise, for the B — > D*ij: decays 



02 

Ol 



(0.42 ± 0.06) e 



i63° 



X 



D*TT 



1 - 






^2^3/2 


D*7r 


\+ 


^1/2 






D*TV 



aBD* 
^0 



.mz 



D* 



(23) 



With the form factors given in various models, we are ready to extract ai and 02 from the 
experimental data. The results are shown in Table III and the parameter 02 falls into the 
range of \a2{D7c) \ ~ 0.35 — 0.60 and \a2{D*TT) \ ~ 0.25 — 0.50. Note that the phases of 02/01, 
59° for the Dn system and 63° for D*n, are slightly different from that given in based 
on the preliminary Belle and CLEO data. We see that although |a2(-Dvr)| and |o2(-D*7r)| 
agree to within one standard deviation, there is a tendency that the former is slightly larger 
than the latter. Hence, nonfactorizable effects could be process dependent, recalling that 
the experimental value for B J/^K is |o2(J/^K)| = 0.26 ± 0.02 0. 

Ideally, the parameters Oi and 02 will be more precisely determined if the topologies 
T, C and S can be individually extracted from experiment. Indeed, this is the case for 
charm decays where T, C and S can be determined from D Kit, D Krj and D Krj' 
decays based on SU(3) flavor symmetry and it is found that \T\ : \C\ : \£\ ~ 1.7 : 1.3 : 



1.0 [18|. Hence, the VT-exchange amplitude that receives short- distance and long-distance 



contributions is not negligible at all in charm decay.[| Unfortunately, one cannot extract those 



tFrom |T| one can deduce that ^02/01 =C/T = (0.73±0.05) exp(fl52°) for D ^ PP decays with- 
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TABLE III. Extraction of the parameters oi and 02 from the measured B — > D^*\ rates 
by assuming a neghgible T^-exchange contribution. Note that a2(-Dvr) and a2{D*TT) should be 
multiphed by a factor of (200 MeY/fo) and (230 MeV//D.)> respectively. 



Model 






a2{DiT)/ai{Dn) 


|ai(D*7r)| 




a2{D*iT)/ai{D*n) 


NRSX 
LF 
MS 
NS 


0.85 ±0.06 
0.84 ±0.06 
0.88 ±0.06 
0.93 ±0.06 


0.40 ± 0.05 
0.44 ± 0.06 
0.47 ± 0.06 
0.56 ±0.07 


(0.47 ± 0.05) exp(j59°) 
(0.53 ± 0.06) exp(j59°) 
(0.53 ± 0.06) exp(j59°) 
(0.60 ± 0.07) exp(j59°) 


0.94 ± 0.04 
0.80 ± 0.03 
0.85 ± 0.03 
0.91 ± 0.03 


0.31 ± 0.04 
0.36 ± 0.05 
0.389 ±0.05 
0.44 ± 0.06 


(0.33 ± 0.04) exp(i63°) 
(0.45 ± 0.06) exp(i63°) 
(0.46 ± 0.06) exp(i63°) 
(0.48 ± 0.06) exp(i63°) 



three quark-diagram amplitudes for B decays since the decay ampHtudes of i? — > D^{ri,ri') 
are proportional to {C + £), while D°7r° is governed by (— C + £) [see Eq. @)]. Therefore, 
the quark-diagram amplitudes C and S cannot be disentangled. Nevertheless, an accurate 
measurement of D^{ri,ri') will enable us to test the importance of W^-exchange in B ^ Dtt 
decays. 

In principle, ai can be determined in a model-independent way from the measurement of 
the ratio of the decay rate of color-allowed modes to the differential semileptonic distribution 
at the appropriate [p!9|: 



S<" . /'^ - = S.'alfHV.,?Yl'\ (24) 

' dB(B° ^ r>W+C-p)/<ig2 , , ' ' " ' 

where Vij is the relevant CKM matrix element and the expression of Y;|*^ can be found in [Q. 
Since the ratio S'^*'* is independent of Vcb and form factors, its experimental measurement can 
be utilized to fix oi in a model-independent manner, provided that vj^'' is also independent 
of form-factor models. Based on the earlier CLEO data, it is found that ai(Z)7r) = 0.93 ±0.10 
and ai{D*TT) = 1.09 ± 0.07 Needless to say, the forthcoming measurements from BaBar, 
Belle and CLEO will enable us to extract the model-independent ai more precisely. Note 
that QCD factorization predicts ai(D^*^'K) ^ 1.05 in the heavy quark limit [pO| . 

Assuming a2{D^*'^ri^'^) = a2{D'^*^'n) we see from Table IV that the predicted branching 
ratios of B^ — > D''*^^t] are consistent with experiment. Note that the predicted rates of 
D^*^^(r],r]') are the same for LF, MS and NS models since a2{D7i) F^^'^ (m'jj) is model inde- 
pendent [see Eq. (^)] and the form factors F^^° and F^^^ are assumed to be proportional 
to F^'" in these models. 



out making any assumption on Vl^-exchange, to be compared with the value (1.05ib0.05)exp(il49°) 
obtained in [^ by neglecting W^-exchange. 
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TABLE IV. Predicted branching ratios (in units of 10 ^) of S D^*^^{ri,ri') in various 
form-factor models by assuming 02(1?''* = a2{D^*'^Tr). 



Decay 


mode 


NRSX 


LF 


MS 


NS 


Experiment [1] 






1.43 ±0.24 


1.69 ±0.28 


1.69 ±0.28 


1.69 ±0.28 


1.4t[]:^±0.3 






1.20 ±0.29 


1.41 ±0.35 


1.41 ±0.35 


1.41 ±0.35 


2.0t[j|±0.4 






0.89 ±0.15 


1.05 ±0.18 


1.05 ±0.18 


1.05 ±0.18 




b"^ 




0.72 ±0.18 


0.85 ±0.21 


0.85 ±0.21 


0.85 ±0.21 





B. Effective parameters and a*^ 

Thus far we have assumed that quark- diagram topologies include all strong-interaction 
effects including FSIs. It is equally well to take a different point of view on the quark-diagram 
topologies, namely, their original forms do not include FSIs from the outset. In this case, 
there is no relative strong phase between the isospin amplitudes A1/2 and given by Eq. 
(p!7|). Next, one puts isospin phase shifts into Eq. (0) to get 

^ /^°7r°)Fsi = ^l\Ai/2e''-" - ^1^3/26^^3/., (25) 
A{B- D\-)fsi = V3Ay2e''-'^\ 

where the subscript "FSI" indicates that the physical amplitudes take into account the effects 
of FSIs. This is motivated by comparing the experimental results with the calculated isospin 
amplitudes under the factorization approximation. Neglecting inelastic scattering, one can 
then extract the coefficients a'f^ from a comparison of the measured and calculated isospin 
amplitudes It is straightforward to show that 

^ D\\si = A(B° ^ D%^) + (e^(^i/2-^3/2) _ 1) ^ 

A(B^ _ D+7r-)Fsi = AiB" ^ D+TT-) ± ^^-^ + 3^ (e^(<5i/2-<53/2) _ 1^ ^ (26) 

where we have dropped the overall phase e^^^/^. The quark-diagram amplitudes T, C, E 



in Eq. (|26| ) have the same expressions as before except that ai 2 in Eqs. (P and (0) are 
replaced by the real parameters a^'l. The latter do not contain FSI effects and are defined 
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for (5i/2 = (^3/2 = fT^.Q In other words, the parameters are defined when FSIs are not 
imposed to the topological quark diagram amphtudes. 

The isospin phase difference in Eq. (p6D is 29° for both B Dir and B D*7c. It is 
easily seen that af^ /af^ is determined from the second line of Eqs. (^2]) and (p3D but without 
a phase for the ratio y4i/2/(v^^3/2)- For example, af^/af^ for B Dn is given by 



,eff 



af 



Dn 



1 - 


Al/2 






D-K 


\+ 


^1/2 




^2^3/2 


Dn 



Dn ml -ml 
fo ml - ml Ff 



Bn 



,m 



(27) 



D) 



The results are shown in Table V. Obviously af^ /af^ and are smaller than the previous 
solution. 

TABLE V. Extraction of the parameters af^ and from the measured B D^*^-it rates. 
Note that af{DTT) and af{D*Tr) should be multiplied by a factor of (200 MeY/fn) and (230 
MeV//z)*), respectively. 



Model 




afiDTT) 


af{DTT)/af{DTT) 


af{D*TT) 


af{D*7T) 


af{D*7r)/af{D*TT) 


NRSX 


0.88 ± 0.06 


0.23 ±0.08 


0.26 ±0.09 


0.97 ± 0.04 


0.16 ±0.04 


0.17 ±0.04 


LF 


0.87 ±0.06 


0.25 ±0.09 


0.29 ±0.10 


0.83 ± 0.03 


0.18 ±0.05 


0.22 ±0.05 


MS 


0.91 ± 0.06 


0.27 ±0.10 


0.30 ±0.10 


0.87 ±0.03 


0.20 ±0.05 


0.23 ±0.06 


NS 


0.96 ± 0.06 


0.32 ±0.12 


0.34 ±0.11 


0.94 ± 0.04 


0.22 ±0.06 


0.24 ±0.06 



C. Comparison 

We are ready to compare the above two different types of approaches. In the type-I solu- 
tion, D^*^^TT^ rates are accommodated because of an enhanced |a2(-D^*''°7r)|. The branching 
ratio of D'^*^^n^ is not over-estimated owing to a relative strong phase between ai and a2. 
In the type-II solution, although af^ is smaller than the magnitude of 02, the Z}(*)°7r° states 



■^The distinction of hard and soft FSI phases in principle cannot be done in a systematical way. 
For example, a sizable "hard" strong-interaction phase for 02 in B ^ vrvr decay is calculable in the 
QCD factorization approach. However, 02 is not computable for B Dn and hence its strong 
phase is most likely soft. 
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gain a feedback from D^*^^'7i~ via FSIs. | More precisely, elastic FSIs will enhance the decay 
rate of D°7r° by a factor of about 3 and suppress D^n^ slightly. 

It has been realized that the isospin analysis proves to be useful only if a few channels are 
open as the case of two-body nonleptonic decays of kaons and hyperons. The isospin phases 
there (or decay amplitude phases) are related to strong-interaction eigenphases (for a recent 
discussion, see [0). For example, one can identify the isospin phase shift in if — > tttt with 
the measured tttt strong- interaction phase at the energy -^i = mj^. However, when there 
are many channels open and some channels coupled, as in D and especially B decays, the 
decay phase is no longer the same as the eigenphase in the S'-matrix. Indeed, the S'-matrix 
in general contains a parameter describing inelasticity. Consider the decay D^7i~ as 

an example. The state D^n^ couples to not only D°7r°, but also D^t], D^t]', Dtttitt channels, 
■ ■ ■, etc. It has been argued that in the heavy quark limit the B decay is dominated by 
multiparticle inelastic rescattering ||23|. As a consequence, even if elastic D^*^7r scattering 
is measured at energies ^/s = tub, the isospin phases appearing in (|I6|) or (^) cannot 
be identified with the measured strong phases. Moreover, the isospin amplitudes are not 
conserved by inelastic FSIs. Therefore, the isospin analysis presented before should be 
regarded as an intermediate step for describing physical decay amplitudes. 

Nevertheless, the isospin decomposition of B ^ Dn amplitudes in Eq. (^) or (|25|) is 
still valid. The isospin analysis is useful in some aspects. First, it provides an indepen- 
dent check on the relative phases among three decay amplitudes. Second, the deviation of 
|^i/2/("\/2^3/2)| from unity measures the degree of departure from the heavy quark limit 
||17|| . Third, the deviation of 02 from characterizes the importance of (soft) FSI contribu- 
tions to the color-suppressed quark diagram, recalling that are defined for the topologies 
without FSIs. This point will be elucidated more below. 

As stressed in [|T^ , the topological quark graphs are meant to have all strong interactions 
included. Hence, they are not Feynman graphs. For example, the genuine VT-exchange 
topology in i? — i> Dti decay consists of not only the short- distance Pi^-exchange diagram 
but also the rescattering graph in which 5^ D^n^ is followed by the strong interaction 
process: {D^7i^)j=i/2 — scalar resonances — D^tt^. Likewise, the process with inelastic 
rescattering from the leading T amplitude into D^n^ via quark exchange has the same 



topology as the color-suppressed tree diagram C [24|. Therefore, color-suppressed tree and 



^Recently, it has been suggested in |21| that quasi-elastic scatterings of D^*^P D^*^P and 
DV — > DV , for example, DP = D'^tt" , D^tt'^ , D^7]q, Df K~ , can explain the enhancement of not 



only L»°7r° but also D'^ T] via inelastic rescattering from the class-I mode B D~^tt 
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ly-exchange topologies receive short- distance and long-distance contributions. 

From Tables III and V we see that Ron = \a2{D7r)/af{D7r)\ ^ 1.75 and Rd*^ = 
|a2(-D*7r)/a2^(D*7r)| ^ 1.95. The corresponding quantities in D Kn decays are Rj^^ ~ 2.0 
and Rj^*^ ~ 1.7, respectively. Therefore, although the relative phase 59° (63°) between 

TT is significantly reduced from the phase 150° between 
D'^ — > and -D° K^^*\^ [|l^], the ratio R does not decrease sizably from charm 



to bottom and, in contrast, it increases for the VP case. It is thus anticipated that in both 
D — > Kvr and B Dn decays, the soft FSI contributions to the color-suppressed topology C 
are dominated by inelastic rescattering .0 Since rj and u are isospin singlets, the conven- 
tional isospin analysis of FSIs is no longer applicable to the final states involving rj or uj. The 
fact that the predicted D^*^^ri rates based on the assumption a2{D^*^ri) = a2{D^*^iT) 

are consistent with experiment (see Table IV) supports the notion that FSIs in B decay are 
indeed highly inelastic. 



IV. DISCUSSION AND CONCLUSION 

Beyond the phenomenological level, it is desirable to have a theoretical estimate of 
a2{D7c). Unfortunately, contrary to the parameter ai(D7r), a2{DTx) is not calculable in 
the QCD factorization approach owing to the presence of infrared divergence caused by the 
gluon exchange between the emitted meson and the {B^ti^) system. In other words, 
the nonfactorizable contribution to 02 is dominated by nonperturbative effects. Neverthe- 
less, a rough estimate of 02 by treating the charmed meson as a light meson while keeping 
its highly asymmetric distribution amplitude yields a2{DTi) 0.25 exp(— ^40°) [Q. Evi- 
dently, large power corrections from long-distance FSI effects are needed to account for the 
discrepancy between theory and experiment for a2{DT[). The rescattering contribution via 
quark exchange, D^ti^ D^tc^, to the topology C in B^ ^ D^n^ has been estimated in 



|]27| using p trajectory Regge exchange. It was found that the additional contribution to 
D°7r° from rescattering is mainly imaginary: a2(-D7r)/a2(-Dvr)„ithout fsis = 1 + 0.61exp(73°). 
This analysis suggests that the rescattering amplitude can bring a large phase to a2(-Dvr) as 
expected. 

In QCD factorization, a2(7r7r) or a2{KTT) is found to be of order 0.20 with a small strong 



*The quark diagram V7-exchange in D ^ PP decays and its phase relative to the topological 



amplitude T are dominated by nearby resonances in the charm mass region |25|, as shown explicitly 
in PI. 
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phase (see e.g. [^]). The fact that the magnitude of a2{DTi) is larger than the short-distance 
one, a2{K7r) or af^lDiT), should not be surprising because the former includes all possible 
FSIs, while the latter is defined without long-distance FSIs. In other words, a2(-D7r) include 
many possible long-distance effects. In the language of isospin analysis, we see from Eq. ( pGf ) 
that 

a2{D7i) = af{D-K) ^-^ [e''^^^'^~^-^'^^ - l) , (28) 

where we have neglected IV-exchange and 

Id - [m^) 

It follows from Eq. (|19D and Table V that a2{DiT)/af{D7T) ^ 1.65 exp(56°). It is worth 
remarking that a2{J/ipK) in B ^ J/ipK decay is calculable in QCD factorization; the 
theoretical result \a2{J/ipK)\ = 0.19l|]:}2 ^ is consistent with the data 0.26 ± 0.02 [§. 
Hence it remains to understand why |a2(-Dvr)| is larger than |a2(-D*7r)| and \a2{J/'4'K)\ or 
why (soft) final-state interaction effects are more important in Dn, D*7i than in J/ipK final 
states. 

To conclude, the recent measurements of the color- suppressed mo des W D(*)°7r° imply 
non-vanishing relative FSI phases among various B Dn decay amplitudes. Depending 
on whether or not FSIs are implemented in the topological quark-diagram amplitudes, two 
solutions for the parameters ai and 02 are extracted from data using various form-factor 
models. It is found that 02 is not universal: |a2(-D7r)| ~ 0.40—0.55 and |a2(-D*7r)| ~ 0.30—0.45 
with a relative phase of order 60° between ai and 02. If FSIs are not included in quark- 
diagram amplitudes from the outset, we have af^i^Dn) ~ 0.23 — 0.32, a'^{D*'n') ~ 0.16 — 0.22 . 
The large value of |a2(-D7r)| compared to a'^i^Dn) or naive expectation implies the importance 
of long-distance FSI contributions to color-suppressed internal VT-emission via final-state 
rescatterings of the color-allowed tree amplitude. 
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